1. IntroductionHow the energy, initially excited in a particular mode, flows into other nonlinear coupled modes in a chain of classical harmonic oscillators (Fermi–Pasta–Ulam model (FPU)[1]) is still not completely understood.[2] The surprising FPU recurrence of energy, which is opposite to the original expectation,[3] actually gives rise to new phenomena, the integrability of nonlinear equations and the dynamical chaos.[2, 4, 5] Even so, people are continuing to look for an energy sharing by modifying the FPU model to higher dimensions or coupling forms.[2] Once the nonlinearity is over some threshold, the FPU model shows strong energy sharing or energy equipartition among normal linear modes.[4, 5] Linearly dependent coupled frequencies help energy sharing even for a very weakly coupled oscillator system.[6] Further investigation shows that the special initially excited mode under resonance conditions,[6] elastic collision interaction[7] or strong (or weak) chaos conditions,[4, 5] leads to strong energy sharing with other coupled normal modes. So far, a statement is reported:[2] energy may flow into other modes with regular behaviors, initially, excited to a higher mode, in some initial period of time (introduction period), and then a stochastic exchange among the modes comes into play, with a practical irreversibility of motion.
Due to the consistency between the Gross–Pitaevskii (GP) equation (governing the basic feature of BECs) and FPU-β model,[8] the investigation of the FPU problem or energy sharing process with BECs has been theoretically suggested almost twenty years ago[9, 10] and its time evolution has been experimentally attempted to be observed.[11] Based on the analytically completed eigenstates for a perfect box of the finite size, corresponding to the normal modes in the FPU problem, one criteria for the appearance of the equilibrium state is reported in Ref. [9]. A transition from a regular dynamics to an irregular one coincides with the onset of the Tonks–Girardeau regime.[10] However, it is surprising that no clear evidence was found in Ref. [11], even though a very strong nonlinear interaction (within Tonks–Girardeau regime) was considered.
In Ref. [11], wherein a relatively long-lived oscillation of each BECs trapped in a harmonic trap was reported in the position and momentum spaces, we suggest that each component BECs (with given momentum
in Ref.[11]) trapped in a harmonic trap could be considered as a normal mode of the FPU problem. With this, the nonlinear couplings between the models are realized by the collision interaction between different components of BECs, and then the dynamics of this system is governed by the well-known coupled GP equations. To understand the energy sharing process in this system, the Husimi distribution function (one quantum distribution function),[15, 16] the Fourier spectrum of time evolution of energy,[6, 17] and the statistical entropy[9, 18] are calculated. Initially, for the excited one of the two components of BECs, its energy flows into the initial static mode with Rabi-like oscillation in the case of the weak nonlinear coupling. With the increase of the nonlinear coupling, the time evolution of energy is modified from quantum beating to chaos dynamics after a critical time. When the nonlinear coupling is strong enough, a strong energy sharing is found by checking the time average of energy, Husimi distribution function, and statistic entropy.
In the following, we present our model, the equations, and the nonlinear coupling interaction in Section 2. In Section 3, we analyze the energy exchanging processes of two-component BECs by investigating the Husimi distribution function, the Fourier spectrum of time evolution of energy, and the statistic entropy. Section 4 is devoted ordinarily to the case of three components. A conclusion is given in Section 5.
2. Model and nonlinear coupling forceA general description of n-component condensate mixture with a one-dimensional (1D) harmonic trap
can be given by the following coupled GP equations within the mean-field theory:
| (1) |
with
and
, where, for simplicity, we assume that all components of atoms have the same mass
m and the number of atoms of each species is fixed and cannot be converted into a different component. Here, the length and the energy have been rescaled in the units of harmonic trap character length
and energy
. The atomic interactions are described by
gij, which are proportional to the two-body
s-wave scattering length
aij, and describe the intra-nonlinear parameters (
i =
j) and inter-component nonlinear coupling (
), respectively. The dynamics of each component in Eq. (
1) performs a harmonic oscillation around its balance position (
x = 0) due to the harmonic trap
V(
x) without the inter-component nonlinear coupling (
). An interesting energy exchanging process is induced by the inter-component nonlinear coupling, which is defined by the scattering lengths between different components and their relative positions
, where
is the average with the
i-th time-dependent component. As shown in Fig.
1, the Gaussian-like interaction energy (defined as
), depends on
, accompanied by an antisymmetric force (
) when
. The repulsive interaction slows down their approaching process, but speeds up their departing process, which is similar with the elastic scattering interaction except for at very short distance.
In this paper, we focus on the energy sharing process, which could be a quantum analogy with the well-known FPU problem studied for more than 50 years.[1, 2] Compared with the initial FPU model and other existing models in BECs,[9, 10] some points should be reminded: i) the same as that in Refs. [1] and [2], the total number of oscillators is an external parameter, but not the same as that in Refs. [9] and [10], wherein it depends on the initial chosen conditions; ii) the intra-nonlinear term is constant, quite different from that in Ref. [1] (frequency-dependence) and Refs. [9] and [10] (dependence on model integral); iii) the modes considered here, are the physically distinguishable components of BECs, but are not the collective motion as before.
3. Rabi, quantum beating and no recurrence energy sharingStarting with the simplest case governed by Eq. (1), the two components (or species), which may be created by imprinting different momenta on one BECs shown in Ref. [11], are trapped in a harmonic trap V(x), and then all of the intra-nonlinear and inter-coupling parameters in Eq. (1) are the same and are denoted as g. Initially, the two components are exactly the same except for
with zero momentum, while
with a nonzero momentum. It is easy to prove that
performs the oscillated behavior when the inter-component nonlinear coupling (
) is turned off, while
keeps static at the bottom of V(x). Turning on the inter-component nonlinear coupling (
), how initially excited energy of
flows to or exchanges with that of
attracts our attention.
Before going into the details, two important quantities are introduced. One is the energy (Ei(t)) for the i-th component, which is defined by
| (2) |
where
and the time-dependent wave function
can be obtained by numerically solving the coupled GP equations (Eq. (
1)) by the second-order splitting-operator method. The other is the Husimi distribution function, which is one kind of quasi-probability distribution function governed by Schrödinger’s equation,
[15, 16] and provides a convenient way to explore the complex dynamics of a quantum system in phase space.
[16] Using the characteristic units of a harmonic trap, Husimi distribution function can be written as
| (3) |
where
is the coherent state in real space,
| (4) |
The same as in Ref. [
16], we have averaged the Husimi distribution over ten-thousand times after the total observation time (
and
in this paper), where
T0 is the period of the oscillation of excited nonlinear component at
.
As shown in Fig. 2, with the increase of the value of g, the initially excited energy flows into the nonlinear coupled component as a Rabi-like oscillation, a quantum beating oscillation, and then finally approaches to a strong energy sharing, depending on how strong the nonlinear coupling (g) is. The complex dynamics of energy evolution E(t) shown in Fig. 2 at g = 9,10 is consistent with the general picture described in Ref. [2]: A beating-like oscillation is observed when
, but one more complex behavior is found when
, where
, which is around
(
) for g = 9 (10) in Fig. 2, is called the introduction time and the inverse with the strength of inter-component nonlinear coupling g.
With the increase of the nonlinear coupling strength from g = 5 to 9, the energy exchanging between two components performs from one clear regular (Rabi-like) behavior to a more complex behavior, as shown in Fig. 2. The corresponding average Husimi distribution functions are shown at the top of Fig. 2. The blue-dotted lines of E(t) show that the total energy of this two-component system is conserved during the observation time. In the weak nonlinear coupling (for example,
), a Rabi-like energy exchanging can be found, and its frequency is proportional to the nonlinear parameter g. The corresponding Husimi distribution function shows that those two components of BECs share the initial excited energy and develop almost the same region in the phase space. Further increasing g, one more frequency emerges, and the performance of the energy exchanging is the same as a beating-like oscillation, while the behavior of Husimi distribution function does not change any more. A critical time (
) is found for g = 9, over which a relatively large amplitude collapses to a relatively small one. This critical value can be very short (
) when g = 10. Meanwhile, the shape of the averaged Husimi distribution function is modified to circles, which may reveal the completion of the energy sharing process.
To thoroughly understand the dynamic complexity of the time evolution of energy shown in Fig. 2, we make Fourier analysis for
and plot the extracted frequencies in Fig. 3. In Ref. [6] and [17], a strong energy sharing (or an equilibrium) could be distinguished by checking its frequency spectrum. The Rabi-like and quantum beating-like oscillations of Fig. 2 are supported by a single or two primary frequencies, shown in Figs. 3(a) and 3(b) for weak nonlinearity. The time evolution of
can also be fitted by the following equation:
| (5) |
where
i denotes how many primary frequencies
in Fig.
3. The
A,
t0, and
B are the fitting parameters. For example, the weak nonlinearity
g = 5, we have only one primary frequency
, and two frequencies
and
for
g = 7. The corresponding fitting behaviors by Eq. (
5) can be found in Figs.
4(a) and
4(b). The two different behaviors can also be read from Fourier spectrums. For
g = 9 and
, we can still fit its time evolution of
(see Fig.
4(c)) by Eq. (
5) with
and
(as the red solid curve shown in Fig.
3(c)). It means that the time evolution of
can be considered as a quantum beating-like one. Meanwhile, in the case of
, we cannot fit
by Eq. (
5) with only two frequencies (see Fig.
4(d)). Actually, we also see a relatively broad range of frequencies around
in Fig.
3(c), and this implies that a strong energy sharing appears. A similar characteristic in frequency spectrum is also found in Fig.
3(d) for
g = 10.
The equilibrium or strong energy sharing can also be proved by checking the statistical entropy S,[2, 9] which is defined as
| (6) |
where
Pi is defined as the probability
, where
is the time average of energy of every component and
is the total energy of all components. The corresponding
Pi of process in Fig.
2 and
S are calculated and listed in Table
1. It is easy to see that with the increase of the nonlinearity,
Pi approaches to 0.5 and the entropy
S to
.
Table 1.
Table 1.
| Table 1.
Probabilities P and statistical entropy S of two components.
. |
4. Energy sharing of three componentsIt is interesting to extend the above investigation to more components, for example, three components.
Similar to that in Section 3, the nonlinearity is increased for involving three components. Now, there are at least six nonlinear parameters gij for
. With the increase of the nonlinearity for all
shown in the first three columns of Fig. 5, Rabi-like and strong energy sharing, except for quantum beating, are found. It is easy to note that the Rabi-like energy exchanging can exist through three components, but with different amount. Since we excite initially the third component (
), but the other two components remain static, we have two equal and small amplifications for
in the case of small nonlinearity. With the increase of the nonlinearity, the amplifications (see probabilities
in Table 2) are increasing. It means that the initial excited energy is going to sharing between the three elements. This can be seen from Fig. 6, where more broad frequency spectrums are displayed with g increasing. To show one equilibrium case, we have to consider one more complex nonlinear coupling, where
,
,
and
. The strong energy sharing can be read from statistical entropy S shown in Table 2, where
.
Table 2.
Table 2.
| Table 2.
Probabilities P and statistical entropy S of three components.
. |
5. ConclusionIn this work, we numerically investigated the energy sharing process in multi-component coupled GP equations with a harmonic trap, with help of the second-order splitting-operator method. Clearly, three different energy sharing behaviors, from regular to chaotic, have been shown for two- and three-component BECs with a harmonic trap. With the increase of the nonlinear interaction, more frequencies, from single to more broad frequency spectrum, are involved into the time evolution of energy. An introduction time is inverse with the nonlinear interaction for the two-component case. After that, a strong energy sharing is determined by checking the Fourier spectrum of the time evolution of energy and also the statistical entropy. The final value of statistical entropy is close to
for two components and
for three components, so we would like to say that our system is in an equilibrium state. Considering the situation in experiment,[11] where they initially created two-component BECs with opposite momenta (
) and realized the dynamics for many times, this actually already is one energy equilibrium state in the point of view of energy. Therefore, even though they could increase the nonlinear interaction to a very strong value (
) and observe a long oscillation of dynamical interaction, the energy sharing evidence cannot be found. Based on this reason, we suggest a similar experiment with only one-component excitation. It is a crucial factor that the two components of the BECs have different initial energy. Based on the above discussion, the experiment conditions of our project are qualified to observe the energy sharing, and an evidence of energy sharing could be experimentally supposed to be obtained in the near future.